Orbital reformation with vanadium trimerization in d 2 triangular lattice LiV0 2 revealed by 51 V 

NMR 
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L1VO2 is a model system of the valence bond solid (VBS) in 3d 2 triangular lattice. The origin of the VBS 
formation has remained controversial. We investigate the microscopic mechanism by elucidating the d orbital 
character via on-site 5I V NMR measurements in a single crystal up to 550 K across a structural transition 
temperature T c , The Knight shift, K, and nuclear quadrupole frequency, Sv, show that the 3d orbital with 
local trigonal symmetry are reconstructed into a d yz d a orbital order below T c , Together with the NMR spectra 
with three-fold rotational symmetry, we confirm a vanadium trimerization with d-d cr bonds. The Knight shift 
extracts the large Van-Vleck orbital susceptibility, x WW = 3.6 x 10~ 4 , in a paramagnetic state above T c , which 
is comparable to the spin susceptibility. The results suggest that orbitally induced Peierls transition in the 
proximity of the frustrated itinerant state is the dominant driving force of the trimerization transition. 

PACS numbers: 71.30.+h, 76.60.-k, 75.50.-y, 75.40.Gb 
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Valence bond ordering in two- and three-dimensional sys- 
tems is a manifestation of geometrical frustration. 1 Across 
Mott transition, a valence bond solid (VBS) phase resides 
nearby unconventional metal 2 - 3 and superconducting phases 3 
in triangular lattice systems. The VBS transition can accom- 
pany by d orbital orders with versatile textures, depending on 
lattice geometry and band filling. 4-9 However, d orbital orders 
governed by frustrated intersite interactions are not apparent 
from lattice distortions and need to be detected by microscopic 
probes. Here we show a tool for detailed characterization of 
d 2 orbital ordering in a model material LiVC>2 with the trian- 
gular lattice. 

ABX2 (A: alkali metal; B: transition metal; X: chalcogen 
ion) with an ordered lock-type structure is a prototype of 
triangular-lattice antiferromagnet. 9 For Ti 3+ and V 3+ , ?2# or- 
bital degrees of freedom play crucial roles in determining 
magnetic ground states. Among ABX2, LiV02 has been ex- 
tensively investigated as a candidate of VBS. 1014 The crystal 
structure is comprised of the triangular lattice of V 3+ (3d 2 ) 
ions in a rhombohedral R3c structure [Fig. 1(a)], A trig- 
onal elongation of VC»6 along the c axis splits the t%g sub- 
level into a lower a\ g singlet and an upper e' g doublet [Figs. 
l(b,c)]. 15 Calculated orbital occupations are almost equiva- 
lent in the a\ g and e' g orbitals due to the small trigonal field 
of 25 meV compared to the bandwidth. 15 Below a first-order 
structural transition temperature T c ~ 450-500 K, a super- 
lattice V3a x -\/3fl x 2c [Fig. 1(a)] appears in x-ray and 
electron diffraction measurements. 11-14 Recent pair distribu- 
tion function analyses of synchrotron data propose lower- 
symmetry lattice distortions depending on the layers^ Al- 
though the way of displacements of vanadium ions is com- 
patible to trimerization j 12 i 13 i 16 the orbital character has not 
been accessible by the diffraction measurements due to Li de- 
ficiency and crystal twinning. 

The magnetic susceptibility of LiVC»2 drops at T c with spin 
gap of several thousands Kelvin estimated from the V nu- 
clear spin-lattice relaxation rate T^ l £U& Hence the ground 
state is most likely nonmagnetic without residual paramag- 



netic spins. The spin-singlet state in the d 2 triangular lattice 
has been theoretically explained as a spin-Peierls transition 
accompanied by a d yz d zx orbital order in a Heisenberg model. 5 
However, the conducting in-plane resistivity, ~ 0.1£2cm, and 
the constant magnetic susceptibility 12,14 above T c may be in- 
compatible to a strongly localized picture. Since the ground 
states are chemically tuned from a metallic LiVSe2 (Ref. |2) 
to an antiferromagnetic insulator NaVC»2 (Ref. Il9l) . a spin- 
singlet phase appear to be located in an intermediate param- 
eter region as a function of electron correlation. A Hub- 
bard model calculation explains such behavior by involving 
a strong trigonal crystal field and also modifies the ground 
state for the Heisenberg limit to be a ferrimagnetic order. 20 
Thus, the origin of the spin-singlet state in OVO2 remains 
controversial and to be elucidated via direct observations of 
the orbital structure. 

In this work, we report detailed 51 V NMR studies on a 
single crystal of LiVC<2, which provide microscopic insights 
into the spin-singlet transition. Beyond previous 51 V and 7 Li 
NMR studies for powder samples, which lack information of 
3d orbitals^^^l we determine d orbital characters by inves- 
tigating anisotropic tensors of Knight shift, K, and nuclear 
quadrupole splitting frequency, 5v. The tensors satisfy an or- 
bital order with the three-fold symmetry in the spin-singlet 
state, which promotes cr-bonding vanadium trimerization. We 
also decompose the orbital and spin Knight shifts, which pro- 
vides the significant residual orbital degeneracy. 

Single crystals of LiVC»2 were obtained from a LiB02-Li20 
flux method in a sealed quartz tube by using LiVC»2 pow- 
ders prepared by a solid state reaction of I^COi and V2O3. 14 
51 V NMR spectra were obtained from the spin-echo measure- 
ments in a constant magnetic field Bo = 6. 105 T and a temper- 
ature range 290-550 K. The angular dependences of 5v and 
K were obtained for the single crystal in rotation planes in- 
cluding the principal axes by utilizing a two-axis goniometer 
at 300 K. High-temperature experiments were performed for 
aligned single crystals with the uniaxial rotation in an oven, 
where the temperature was monitored by a thermocouple near 
the sample. 
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FIG. 1 : (Color online) (a) Layered triangular lattice of LiV02 viewed 
from the c axis with an orbital texture obtained from the present 
study, where only two of four leaves in d yz and d zx orbitals are shown 
for clarity. Thick solid and dotted lines represent unit cells of an 
original R3m lattice above T c and a y/3a x ^3a superlattice below T c , 
respectively (Refs[T33). (b) VC>6 unit trigonally elongated along c 
axis above T c . (c) d levels in a trigonal field: lowest a\ g singlet, mid- 
dle e' g doublet, and higher e g doublet, (d) Minority d xy orbitals that 
dominates the anisotropy of NMR spectra below T c . The x, y, and z 
directions are taken parallel to the VO bond directions and identical 
to Y - X, X + Y, and Z axes, where X, Y, and Z are the principal 
axes of the local symmetry that determines the K and Sv tensor. The 
b" axis is defined perpendicular to the a and c axes. The b' axis is 
located at 30° from the a axis. 



NMR is useful as a local probe of d and p orbital occu- 
pations in transition metal 22-25 and organic compounds^ via 
anisotropic electric and magnetic hyperfine interactions. For 
tz g orbitals weakly bonded to ligand p orbitals, the symmetry 
axis of K and Sv reflects occupied 3d orbitals. The method can 
be applied in the spin-singlet state. 27 For 5I V nuclear spins, 
/ = 7/2, electric quadrupole interactions under the electric 
field gradient (EFG) split NMR spectra with an interval fre- 
quency 6v. 6v consists of the on-site d orbital and the outer 
ligand terms: Sv = 5v° n + 6v° m , where the components of 5v° n 
are expressed by 22 



3e 2 Q 
2/(2/- 1) 



(1) 



where e is the electron charge, Q the nuclear quadrupole mo- 
ment, -0.05 x 10~ 24 cm 2 for 51 V, and (r 3 ) the V 3+ radial ex- 
pectation value with a reduction factor of 0.8 compared to the 
free ion value, £ = l/21,and#y = j(LjLj+LjLj)-^L(L+l)5ij 
(L: orbital angular momentum). 6v°J is proportional to the or- 
bital polarization reflected in qij. For instance, a full polariza- 
tion of dxy gives the principal components q xx = q yy = — 1 and 



q zz = 2. 22 - 26 On the other hand, 



<5v out is given as 22 
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shell polarization, and Vy the EFG tensor of surrounding 
ionsi^ 2 . Magnetic hyperfine interactions between 51 V nuclear 
spin and electrons produce Knight shift, where the compo- 
nents are derived from a spin-Hamiltonian expressed in a 
second-order approximation as^ 



Kij = P(kS u + 3&/U + 2AA i7 )^ + 2<Pxlf 



(3) 



by using f = 2fiB~y n ti(r 3 ) the Bohr magneton, y n : the 
nuclear gyromagnetic ratio), a dimensionless parameter k, a 

... . i a s-. <0|L/|n)(n|L,|0> 

spin-orbit coupling parameter A, Ay = 2j K , 



E„-E 



where is the Sternheimer antishielding factor due to core- 



(the 

ground state energy: Eq, the excited one: E„), the spin suscep- 
tibility^., and the Van-Vleck susceptibility = 2Afy/ 2 Ay. 
The term proportional to x]j m Eq. (3) consists of isotropic 
core polarization, 3d dipole interactions, and a cross term with 
spin-orbit coupling. 

To determine Kjj and dv/j in the spin-singlet state, 51 V NMR 
spectra were measured by rotating the single crystal around 
the a, b* , and b axes [defined in Fig. 1(d)] at 300 K, as shown 
in Fig. 2. A single line at Bq\\c indicates that all V sites are 
equivalent, simultaneously satisfying 5v = vzzO c °s 2 #o~l) = 
(6>o = 54.7°: the magic angle). The spectra split into sev- 
eral lines at the other angles due to appearance of inequivalent 
V sites with clear quadrupole splittings: these are assigned to 
two or three sets of seven lines, as traced with fitting curves. It 
indicates symmetry lowering from the original R3m structure 
where all V sites are equivalent at arbitrary Z?o directions. In 
the R3m structure with the local trigonal distortion, the prin- 
cipal Z axis having a maximum 5v, should be parallel to the 
c axis, as known in V2O3. 29 At 300 K, however, we observe 
a 6v maximum around Qq rotated from c to b* axis and ±60 
from a to c axis, where 9o matches the V-O bond direction 
in the VC>6 octahedron. In addition, identical spectral pro- 
files between the cb* and cb' rotations (b' is located at 120° 
inclined from b*) indicate that the observed three V sites are 
connected by three-fold rotational symmetry expected from 
the in-plane superlattice y/3a x y]3a of the x-ray and electron 
diffractions ili— 

We define 5v as an average of the splitting frequency in- 
dependently fitted to the sinusoidal function and plot as a 
function of the rotation angle in Fig. 3. 6v obtained for 
the three V sites are well fitted to the Volkov's formula 28 
and give the principal components, (|<5vxxl, l^vyy!, I^vzzl) = 
(0.15, 0.33, 0.48)MHz, where X, Y, and Z are parallel to x -y, 
x + y, and z in Fig. 1(d) and satisfy the three-fold rotational 
symmetry. The absolute values are good agreement with those 
for the powder sample, (0.15, 0.33, 0.47)MHz. 21 The symme- 
try Z axis parallel to the VO bond z direction in Fig . 1 (d) in- 
dicates that the predominant d (electron or hole) orbital lies 
down perpendicular to the Z axis, i. e., has a ^-f 2 (d^) 
character. In the d 2 case, a minority orbital (a hole orbital) 
can contribute to the EFG. Namely, electrons are mainly oc- 
cupied in dyz (d yz ) and dzx (d zx ), and the sign of £v„ should 
be selected as (6v xx , 5v YY , Sv zz ) = (-0.15, -0.33, 0.48)MHz. 
In a point-charge approximation, <Sv out is obtained as 
(6v X x, Sv YY , Sv zz ) = (-0.01, 0.04, 0.05)MHz for y M , based on 
the crystal structure, 16 where the Z axis is close to the c axis 
due to the residual trigonal distortion. The deviation of the 
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FIG. 2: (Color online) Frequency-swept 51 V NMR spectra in LiV0 2 (300 K, B = 6.105 T). Magnetic field is rotated around the (a) a, (b) b", 
and (c) b axes. Dotted guide lines are sinusoidal fitting 28 at the bottom of the peak position by taking the vertical axis as a rotation angle. Blue, 
red, and green lines are assigned to VI, V2, and V3 in Fig. 1(c), respectively. 



Z axis from the experimental results suggests the negligible 
6v° ut contribution likely due to the small 1 - y x in the present 
case. From Eq. (1), the orbital occupation is evaluated as 
dxi-vidxy) : d Y z(d yz ) : d Z x(d zx ) = 0.34 : 0.77 : 0.89. The 
result does not contradict to the d yz d zx orbital order accompa- 
nied with the trimer-type valence bond order proposed theo- 
retically by Pen et al. 5 but shows a significant contribution of 

dxy. 

Since x s already goes to zero at 300 K, the anisotropy in K 
is governed by the Van-Vleck term of Eq. (3). The symmetry 
of the 3d orbital state is reflected in the anisotropic Van-Vleck 
term via the second-order mixing of the ground and excited 
states under the external magnetic field. K defined as a relative 
shift of the central line is also plotted in Fig. 3. The angular 
dependences of K exhibit similar profiles to those of 5v for the 
three rotations. The principal components (Kxx, Kyy,Kzz) = 
(0.46, 0.47, 0.12) ± 0.02% give an isotropic part K iso = (K xx + 
K YY + K zz )/3 = 0.35% corresponding to x™ = 4.9 x 10~ 5 
emu/V-mol, using the orbital hyperfine coupling -2P = 40 
T (Ref. I22I) for V 3+ . The Z axis is identical to that of 6v, 
consistent with the predominant d xv hole state [Fig. 1(c)], 

To gain further insights into the phase transition, 51 V NMR 
spectra were measured across T c at the c axis and 0q, as shown 
in Fig. 4. A low-frequency (67.5 MHz at 300 K) peak with 
a negative shift comes from an impurity phase of L1V2O4. 30 
For Z?olk, a sharp peak is observed at 68.5 MHz and suddenly 
disappears due to fastening of the spin-echo decay time T2 on 
heating above T c ~ 500 K. Such abrupt fastening of T2 except 
for the magic angle is also encountered in correlated param- 
agnetic metals V2O3 and LiV204i^£ Alternatively, multiple 
peaks observed for fioll^o become single above 450 K. A dif- 
ference of T c in the two measurements arises from thermal 
hysteresis of the first-order transition. Such a switching of the 
magic angle indicates a change in the symmetry axis of the 
d orbital across T c . The magic angle of 6q above T c agrees 



with the trigonal symmetry*** 3 * as expected from the crystal 
structure. 

From the resonance line at the magic angles in Fig. 4, we 
obtained the Knight shift of the paramagnetic state as K\ so = 
0.67%, which is larger than that of the nonmagnetic state, 
0.35%. Here K iso is expressed as P(k + 2AK ii )x) i + 1T>xV 
from Eq. (3). Since k + 2AAu would be negative for vanadates 
owing to the predominant core polarization, an increase in K 
above T c indicates the appreciable spin-orbit coupling term 
and the Van-Vleck term 2 l PxJj V ■ The hyperfine cou- 
pling constant (-6 T//zb) 21 is indeed much larger than that of 
V2O3 (-13 T/jUb) 29 - The enhanced;^ Y v comes from the resid- 
ual degeneracy in the three-fold tz g levels, as expected from 
the calculated orbital occupations in the paramagnetic state. 15 
Using <PU + 2AKu) = -6 T^u B (Ref. 17.) and -29 = 40 T/^ B 
(Ref. |H, we obtained K vv = 2.7% and x VV = 3.6 x 10~ 4 
emu/mol, which corresponds to 47% of the bulk magnetic sus- 
ceptibility (7.6 X 10" 4 emu/mol)^ By subtracting x vv from 
the bulk susceptibility, we obtained^ ~ 4.0 x 10~ 4 emu/mol. 
This value is even lower than that of the metallic phase V2O3 
(Ref. |H) and suggests the itinerant nature of the paramagnetic 
phase of L1VO2. 

The above results have revealed the local orbital structures 
in the triangular lattice LiVC>2 across the structural transition 
by the investigation of the electrostatic and magnetic hyper- 
fine couplings. In the paramagnetic state, the orbital state is 
governed by the trigonal crystal field that remains degeneracy 
in e' r The moderate value of the spin susceptibility suggests 
the system close to a metal-insulator crossover or a metallic 
region. 2 Below the first-order structure transition temperature, 
we observed the d xy hole orbital order that makes occupied 
d yz d zx spins inert due to the strong singlet formation via the 
direct d-d bonding. This trend is promoted by the large dis- 
placement of V atoms. The V-V distance (2.56A) is shorter 
than that of V metals (2.61A) and V0 2 (2.65A), 7 suggest- 
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observed residual Cf^Y = 4.9 x 10 5 emu/V-mol) suggests 
the low-lying inter-orbital excitation on the single ion. The 
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FIG. 3: (Color online) Angle dependence of (a,b,c) nuclear 
quadrupole splitting frequency Sv and (d,e,f) Knight shift A* obtained 
for the b*c, ca, and b'c rotations. X and Z are defined by a maximum 
and minimum of K and Sv, respectively. Experimental data are well 
fitted into the Volkov's sinusoidal formula, y = yfj + ~f x cos 2(9 - <jf) 
(y = Sv, K), with fitting parameters yJJ, yj* and for the rotation axes 
a = a, b* , and bM 



ing the metallic bonding formation. Thus, the present orbital 
oder appears in order to gain the intersite d-d exchange en- 
ergy or form the trimer orbital, instead of the Jahn-Teller dis- 
tortion. Unlike the d l spin-dimer insulators such as MgTi2C»4 
and VO2, the orbital polarization in d 2 system indicates that 
each vanadium site has a spin S ~ 1, similar to the Mott 
insulator. 5 However, the spin gap (1500-3400 K)^! in Li VO2 
is too large for the spin-Peierls Mott state where the spin gap 
should be much smaller than the charge gap, but compara- 
ble to the charge gap. 14 Taking these experimental data into 
account, one can regard the system as an orbitally polarized 
band insulator. Namely, the structural phase transition is the 
orbital-assisted Peierls transition. 7 The appearance of the cor- 
related band insulator agrees with the numerical calculations 
of the Hubbard model. 20-31 Although the orbital occupation is 
much larger than the theoretical prediction, 20 the V displace- 
ments that was not taking into account in the calculation may 
enhance orbital polarization. 

In the band insulator without orbital degrees of freedom, 
one expects that a band gap or a bonding-antibonding orbital 
band transition is equivalent to a spin gap. For the band insula- 
tor with the gap of thousands Kelvin, the orbital susceptibility 
should be vanishingly small (^ vv < 10~ 6 emu/V-mol). The 
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FIG. 4: (Color online) Temperature dependence of the 51 V NMR 
spectra measured at (a) the c axis and (b) 54° from the c axis in 
LiVCK. The spectrum pointed by the arrow comes from the impurity 
LiV 2 4 . 



characteristic in the orbitally polarized Peierls state and pro- 
vide important insights into the Mott versus Peierls arguments 
in transition metal compounds such as V02^ 

In conclusion, we addressed the microscopic origin for the 
anomalous spin-singlet formation in LiV02 with the triangu- 
lar lattice. Our orbital-resolved NMR study revealed the d 
orbital occupation across the structural phase transition. We 
observed the d xy hole order with the three-fold symmetry in 
the Knight shift and nuclear quadrupole coupling for the spin- 
singlet state. The rotation of the local symmetry axis gives 
microscopic evidence for the orbital reformation across T c . 
The origin of the large spin gap arises from the direct d-d cr 
bonding between d yz and d zx orbitals in the vanadium trimer. 
In this respect, the phase transition is driven by the strong 
electron-phonon coupling and orbital ordering to relieve the 
geometrical frustration and effectively gain the intersite ex- 
change energy, which can be regarded as the two-dimensional 
Peierls transition. The residual low-energy orbital excitation 
may feature the orbitally-polarized band insulator with mod- 
erate electron correlations. 
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